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ABSTRACT 

Small hydrocarbons, such as C 2 H, C 3 H and C 3 H 2 are more abundant in photo-dissociation regions 
(PDRs) than expected based on gas-phase chemical models. To explore the hydrocarbon chemistry 
further, we observed a key intermediate species, the hydrocarbon ion LC 3 H+, in the Horsehead PDR 
with the Plateau de Bure Interferometer at high-angular resolution ( 6 "). We compare with previous 
observations of C 2 H and 0 C 3 H 2 at similar angular resolution and new gas-phase chemical model 
predictions to constrain the dominant formation mechanisms of small hydrocarbons in low-UV flux 
PDRs. We find that, at the peak of the HCO emission (PDR position), the measured /-C 3 H + , C 2 H 
and c- C 3 H 2 abundances are consistent with current gas-phase model predictions. However, in the 
first PDR layers, at the 7.7 /im PAH band emission peak, which are more exposed to the radiation 
field and where the density is lower, the C 2 H and 0 C 3 H 2 abundances are underestimated by an order 
of magnitude. At this position, the /-C 3 H + abundance is also underpredicted by the model but only 
by a factor of a few. In addition, contrary to the model predictions, /-C 3 H + peaks further out in 
the PDR than the other hydrocarbons, C 2 H and 0 C 3 H 2 . This cannot be explained by an excitation 
effect. Current gas-phase photochemical models thus cannot explain the observed abundances of 
hydrocarbons, in particular in the first PDR layers. Our observations are consistent with a top- 
down hydrocarbon chemistry, in which large polyatomic molecules or small carbonaceous grains are 
photo-destroyed into smaller hydrocarbon molecules/precursors. 

Subject headings: astrochemistry molecular data molecular processes ISM: abundances ISM: 
molecules photon-dominated region (PDR) 


1. INTRODUCTION 

Simple hydrocarbon molecules, such as C 2 H, C 3 H, 
C 3 H 2 and C 4 H are ubiquitous in the interstellar medium 
(ISM). They are easily obse rved in a large var iety 
of sources, from diffuse (e.g. [Lu cas & Liszt 200 (1) to 
dark clouds (e.g. IWootten et all 119801 : ICox et ahl 119891 : 
iMangum fe Woottenlll990[ ). In Photon Dominated Re¬ 
gions (PDRs) they have been found to be almost as 
abundant as in dark, well-shielded clouds, despite the 
strong UV radiation co mpared to the mean interstel¬ 
lar radiation field (e.g.. iFosse et a l. 2000; iFuente et al.l 
l2003t iTevssier et al.l l2004 iPetv et aTl l2005h . In con¬ 
trast to high UV-fl ux PD Rs (y ~ 10 4 — 10 5 relative 
to the Draine field; iDrainel fl978l ). like Mon R2 or the 
Orion Bar, where the observed hydrocarbon abundances 
can roughly be explained with pure gas-phase chemistry 
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models (jCuadrado et al.l l2014f h in low-UV flux PDRs 
(x < 100 ) the current pure gas-phase ch emical models 
fail to reproduce their hi gh ab undances. ITevssier et al.l 
(12001 and IPetv et ahl (j2QQ5h proposed that another 
mechanism producing carbon chains must exist in ad¬ 
dition to gas-phase chemistry. One such possibility, sug¬ 
gested by some laboratory experiments and theoretical 
calculations, is the fragmentation of polycyclic aromatic 
hydrocarbons (PAH) or very small carbo naceous grains 
(VSGs) due to the far-UV ra diation field (jLe Page et al.l 
I2OO3I : iMontillaud et al.ll2013l . and references therein). In¬ 
deed, recent laboratory experiments have shown that 
the far-UV irradiation of interstellar hydrogenated amor¬ 
phous carbon analogs leads to the effi cient pr o ductio n 
of small hydrocarbons, such as CH 4 (jAlata et al.l 120141 ). 
The good spatial correlation between the mid-IR emis¬ 
sion due to PAHs and the distribution of carbon chains 
found in the Horsehead nebula, a low-UV flux PDR 
(x ~ 60), provides some support for a hy drocarbon 
produ ction mechanism starting from PAHs (jPetv et al.l 
l2005h . 

Using single-dish observations, IPetv et all d 2 f)T^ re¬ 
ported the first detection of /-C 3 H + , a key intermediate 
species in the gas-phase formation of small hydrocarbons, 
toward the Horsehead PDR. /-C 3 H + wa s later detecte d 
towa rd the Sgr B 2 (N) molec ular cloud (jMcGuire et al.l 
120131 ). iBriinken et al.l (j2014h measured the millimeter 
rotational spectrum of /-C 3 H + in the laboratory and un¬ 
ambiguously confirmed the assignment of the observed 
lines in th e Horsehead to th e hydr ocarbon ion, LC 3 HV 
Moreover. iBotschwina et ahl (1201 4) performed highly ac¬ 
curate quantum chemical calculations to study this ion 
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TABLE 1 

Observation parameters for the maps. Their projection center is <32000 = 05 /l 40 m 54.27 s , <$2000 — — 02°28 / 00 // . 


Line 

Frequency 

GHz 

E u /k 

K 

Instrument 

Beam 

arcsec 

PA 

0 

Vel. Resol. 
km s -1 

Int. Time 6 
hours 

T s ys 

k (TO 

Noise 

K (Tint.) 

C-C 3 H 2 22,1 — lo,l 

85.339 

6.4 

PdBI/C&D 

6.1 x 4.7 

36 

0.2 

12.0 

- 

0.30 

C 2 H N = 1 - 0° 

87.317 

4.2 

PdBI/C&D 

7.2 x 5.0 

54 

0.2 

12.0 

- 

0.60 

Z-C 3 H+ J = 5-4 

112.446 

16.2 

PdBI/C&D 

6.2 x 5.5 

25 

0.2 

17.3 

150 

0.07 

DCO+ J = 2-1 

144.077 

10.4 

30m/CD150 

18.0 x 18.0 

0 

0.1 

1.5 

230 

0.10 


a Transition: N=l—0,J = 3/2 — 1/2,F = 2 — 1. b On-source integration time scaled to a 6 antenna array. 
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Fig. 1.— Integrated intensity maps of the small hydrocarbons Z-C3H"*", C2H and C-C3H2 lines (upper row), as well as that of the DCO + 
line, the 7.7 fim PAH emission and the H 2 2.12 /xm ro-vibrational line (bottom row). Maps have been rotated by 14 counter-clockwise 
around the projection center, located at ( Sx,5y ) = ( 20 // , 0 // ), to bring the illuminating star direction in the horizontal direction and the 
horizontal zero has been set at the PDR edge, delineated by the red vertical line. The blue and red crosses show the Core and PDR 
positions, respectively. The two gray horizontal lines on top of the /-C3H+ map display the region over which the spectra are averaged in 
Fig. [2] All maps have been integrated between 10.1 and 11.1 kms -1 . 


and found an excellent agreement with the spectroscopic 
constants derived fr om the observations i n the Horse- 
head. More recently. ICnadrado et al.l (|2014l) detected the 
/-CsH + lines up to J = 13 in the Orion Bar, what allowed 
the m to further refine t he rotational constants. In addi¬ 
tion, iMladenoviel ([20141 ) revised the spectroscopic param¬ 
eters of /-CsH + by means of numerically exact rovibra- 
tional calculations. The detection of the small hydrocar¬ 
bon /-CsH + has thus been confirmed, both theoretically 
and experimentally. 

Chemically, CsH + is a gas-phase precursor of the small 
hydrocarbons C 3 H and C 3 H 2 . Observations of the dis¬ 
tribution of C 3 H+ with respect to these hydrocarbons 
could thus be used to constrain the dominant formation 
mechanism of small hydrocarbons. The Horsehead Neb¬ 
ula provides an ideal test-bed because it is viewed almost 
edge-on, providing easy access to the warm surface layer 
of a cloud where C + and therefore C 3 H + is predicted to 
be the most abundant. A s this warm photo-a ctive layer 
is spatially narrow SA lGuzman et al.ll 2012 [ ). high an¬ 
gular resolution is needed to resolve the steep gradients in 


this region. To explore the relationship of CsH + with its 
environment and with neutral carbon chains, we present 
in this letter the first spatially resolved observations of 
Z-C 3 H+. 

2. OBSERVATIONS 

We used the Plateau de Bure Interferometer (PdBI) 
to obtain a 6" angular resolution map of the emission of 
the /-C 3 H + (5 — 4) line at 112.446 GHz. The observation 
parameters are summarized in Table [lj The observa¬ 
tions were carried out in December 2012 , April, May, 
September and October 2013 with six antennas in the C 
and D configurations (baseline lengths between 24 and 
176 m). We observed a 9-field mosaic and used about 37 
hours of telescope time, which correspond to 17.3 hours 
of on-source time scaled to a six antenna array after fil¬ 
tering out low-quality visibilities. The /-C 3 H + line was 
covered with a correlator window of 20 MHz bandwidth 
and 39 kHz channel spacing. The typical precipitable 
water vapor amounted to 6 mm and the typical system 
temperature was 150 K. 
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Spatially resolved /-CsH + emission in the Horsehead PDR 


The PdBI data was calibrated with the GILDAS FV CLIC 
software. The bright quasars 3C84, 2200+420 and 3C279 
were used to calibrate the radio-frequency bandpass, and 
two nearby quasars (0420—014 and 0528+134) were reg¬ 
ularly observed to calibrate phase and amplitude tem¬ 
poral variations. MWC349 was used to derive the abso¬ 
lute flux scale. In order to recover the extended emis¬ 
sion that is filtered out by the PdBI, we observed the 
same region with the IRAM-30m telescope during ^15 
hours of average summer weather in July and October 
2013. The GILDAS/CLASS software was used to process 
the IRAM-30m data and produce the single-dish map, 
which was then combined with the PdBI observations 
in GILDAS/MAPPING, in the same way as described in 
iGuzman et al ] (Imah . 

3. RESULTS 

3.1. Spatial distribution 

Fig. [T| displays the integrated intensity maps, at similar 
angular resolution (6"), of the /-CsH + J = 5 — 4, C 2 H 
N = 1 -0 , J = 3/ 2 — 1/2, F = 2-1 and c-C 3 H 2 2 2 ,i -l 0 ,i 
lines (|Petv et al.ll2005| ). Also displayed are the 7.7 /am 
emission arising from PAHs imaged with ISO at 6" an¬ 
gular resolution (jAbergel et al.ll2003l ) , and the integrated 
intensity maps of the DCO + J = 2 — 1 line o bserved with 
the IRAM-30m telescope (|Petv et all 1200 7) and the H 2 
v = 1 — 0 5(1) li ne obtained with S OFI/NTT at 1" an¬ 
gular resolution (jHabart et al.ll2005l h 

The emission of the C 2 H and C-C 3 II 2 lines is struc¬ 
tured into two successive filaments parallel to the dis¬ 
sociation front. The first one mainly coincides with the 
PDR, where intense emission is seen both in the H 2 ro- 
vibrational line and in the PAH mid-infrared band. The 
second one is located at the position of the cold dense 
UV-shielded core as indicated by the bright DCO + emis¬ 
sion, with a difference between the two species: At the 
peak of the DCO + emission there is a clear, localized 
deficit of the C 2 H emission compared to that of c- CsH 2 . 
In contrast, the /-CsH + J = 5 — 4 line only emits toward 
the UV-illuminated edge, at the surface of the cloud. 
Moreover, the /-CsH + emission reaches the red line in 
Fig-El which traces the edge of the PDR, while the emis¬ 
sion of the two other hydrocarbons is shifted left of this 
line. We note that the observed /-CsH + line is weak (6a 
at the PDR), and hence the clumpy structure is most 
likely an artifact caused by the low signal-to-noise ratio. 
Indeed, due to the low inclination of the source the dirty 
beam has side-lobes which produce some uncertainty in 
the deconvolution of the data. 

Line spectra along the direction of the illuminating 
star, centered at dy = 0", are shown in the upper panel of 
Fig. [21 To increase the signal-to-noise ratios we have av¬ 
eraged the spectra over 20" in the ^-direction, z.e., the 
region between the two gray lines in Fig.Q] In the bottom 
panel the integrated line intensity profiles of the hydro¬ 
carbons are shown, as well as that of the H 2 line and 
the PAH emission. The dashed vertical lines mark three 
characteristic positions in the Horsehead: the Core, cor¬ 
responding to the peak of the DCO + line emission (RA 
= 5 h 40 m 55.61 s , Dec = -2°27'38", J2000), and charac¬ 
teristic of the cold UV-shielded gas (marked by the blue 

10 See http://www.iram.fr/IRAMFR/GILDAS for more informa¬ 
tion about the GILDAS softwares flPetvl 120050 . 


CORE PDR PAH 
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Fig. 2 .— Left: Spectra along the direction of the exciting star 
at the PDR position (dy = 0 in Fig. [T}, averaged over 20" in the 
dy-direction. Right: Integrated line intensities along the direction 
of the exciting star at the PDR position. 

cross in Fig [Tj) ; the PDR, corresponding to the peak of 
the HCO line emission e mission (RA = 5 h 40 m 53.936 s , 
Dec = — 2°28 / 00 // , J 2 QQQ: iGerin et ahll2009[) and charac¬ 
teristic of the warmer UV-illuminated gas (marked by the 
red cross in Fig. [Tj) ; and a position closer to the edge of 
the cloud, corresponding to the peak of the 7.7 /am PAH 
emission (RA = 5 h 40 m 53.6 s , Dec= -2°28 / 1.9 // , J2000). 
We observe a shift between the peak of /-CsH + emission 
compared to that of the other hydrocarbons. /-CsH + 
peaks at Sx ~ 10 ", while the other hydrocarbons peak 
at Sx rsj 15". Although the peak position of the /-CsH + 
line has a higher uncertainty given the lower signal-to- 
noise ratio, the emission profile is clearly broader and it 
extends further out in the PDR compared to the C 2 H 
and 0 C 3 H. Moreover, the /-CsH + integrated line inten¬ 
sity profile correlates better to that of the H 2 line and 
the PAHs emission compared to neutral hydrocarbons. 

3.2. Abundances 

To compare with chemical models, we have computed 
the column densities of C 2 H, C 3 H, C 3 H 2 (both linear and 
cyclic species) and /-C 3 H + at three different positions, 
namely the Core, PDR and PAH positions. These posi- 
tions are slightly different from those used by iPetv et HI 
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TABLE 2 

Column densities and abundances with respect to H nuclei. 



CORE 

PDR 

PAH 

Nu 

6.4 x 10 22 cm -3 

3.8 x 10 22 cm- 3 

6.4 x 10 21 cm -3 

N( C 2 H) 

< 8.8 x 10 13 

(1.3 - 1.6) x 10 14 

(1.5 -5.0) x 10 14 

N(c- C 3 H) 

(0.8 - 2.3) x 10 12 

(2.4 - 7.2) x 10 12 

- 

N(l- C 3 H) 

(1.4 - 4.2) x 10 11 

(0.6 - 1.8) x 10 12 

- 

JV(o-c-C 3 H 2 ) 

(2.3 - 3.0) x 10 12 

(3.8 - 6.9) x 10 12 

(1.2 - 2.7) x 10 13 

Af(p-c-C 3 H 2 ) 

(0.6 - 1.1) x 10 12 

(1.2 - 1.9) x 10 12 

- 

N(l- C 3 H 2 ) 

(0.9 - 2.6) x 10 11 

(0.5 - 1.5) x 10 12 

- 

iV(2-C 3 H + ) 

< 6.5 x 10 10 

(1.6 - 4.8) x 10 11 

(1.4 -4.1) x 10 11 

X(C 2 H) 

< 1.8 x 10“ 9 

(1.9 - 5.9) x 10" 9 

(1.3 - 9.0) x 10" 8 

X(C 3 H) 

(0.9 - 5.1) x 10“ n 

(0.5 - 2.7) x 10“ 10 

- 

X(C 3 H 2 ) 

(2.8 -8.8) x 10“ n 

(0.9 -3.2) x 10“ 10 

(1.5 - 6.6) x 10“ 9 

X(2-C 3 H+) 

< 1.7 x 10“ 12 

(0.2 - 1.4) x 10“ n 

(1.2 - 7.2) x 10 -11 


d2nn^ and iPetv et all ([2012D , and were chosen to take 
advantage of Horsehead WHISPER line survey. For the 
first two positions, we used the single-dish deep inte¬ 
grations obtained in the line survey and included beam 
dilution factors obtained from the higher-angular reso¬ 
lution PdBI observations when available. As no line 
survey has been made at the PAH position, we de¬ 
rived the abundances directly from the PdBI observa¬ 
tions. For C 3 H, we assumed the same spatial distribu¬ 
tion of c- C 3 H 2 . We used the non-LTE radiative transfer 
code RADEX (|van der Tak et aI1l2007[ ) for those species 
with known collisional rate s, z.e., C 2 H and C -C 3 H 2 (or¬ 
tho and para), taken fr om iSpielfiedel et al.l (J2012) and 
I Chandra &; Kegel] (j2QQQh . respectively. The gas density 
was fixed to 6 x 10 4 cm " 3 in the PDR, 10 5 cm " 3 in 
the dense core (jGerin et al.l 120091 : IPetv et ahll2QQ7f ) and 
(5 — 10) x 10 3 cm -3 in the PAH position. The kinetic 
temperature was left as a free parameter, the best fits 
being consistent with previous estimates of ~ 60 K and 
~ 20 K in the PDR and dense core, respectively. For 
C 3 H and /-C 3 H 2 we constructed rotational diagrams. We 
note that C 2 H and /-C 3 H + are not detected at the core 
position in the PdBI maps. We thus consider their de¬ 
rived abundances as upper limits. The inferred column 
densities and abundances with respect to total hydro¬ 
gen nuclei are summarized in Table 2. The errors in the 
abundances take into account a 50% uncertainty in the 
assumed TVh, which are inferred from the 1.2 mm dust 
c ontinuum emission . 

IPetv et aTI (| 2012 l) derived the /-C 3 H + abundance at 
the PDR position using single-dish observations at ~ 25" 
angular resolution, and computing a beam filling factors 
by assuming that the /-CsH + emission filled a Gaussian 
filament of ~ 12 " width in the Sx direction. This as¬ 
sumption was based on the morphology of the HCO line 
emission, which shows a filament of roughly this width. 
The new PdBI observations at 6 " angular resolution show 
that the /-C 3 H + emission indeed arises from a ^ 12 " fil¬ 
ament, where the /-C 3 H + J = 5 — 4 line is ^ 3 times 
brighter than what was observed with the 30m. 

3.3. Chemistry 

In order to test our current knowledge of the gas-phase 
chemistry of hydrocarbons, we have used an updated ver¬ 
sion of the one-dimensional, ste ady-state photochemical 
code from ILe Petit et all 00061) . The same model was 


used in IPetv et al.l ([ 201 2 b except we have now introduced 
the formatio n and destruction (by H 2 ) rates of CsH + 
measured bv iSavic fe Gerlichl ( 2005 ). which have an in¬ 
verse temperature dependence. For the physical condi¬ 
tions in the Horsehead (nn ~6xl0 4 cm- 3 ,T kin -60 K), 
this results in higher abundances of C 3 H and C 3 H 2 . The 
model includes grain surface reactions for the formation 
of H 2 and other species, such as H 2 CO and CH 3 OH, but 
only gas- phase reactions f or the formation of hydrocar¬ 
bons (see IPetv et al.lI 2012 L for a detailed description of 
the model). 

In high UV-flux PDRs (T kin ~ 100 — 500 K), the forma¬ 
tion of hydrocarbons starts with the formation of CH + 
([Cuadrado et al.ll2014f) thro ugh the very endothe rmic re¬ 
action C + +H 2 -a CH + + H ([Agundez et al.ll2Q10h . In the 
Horsehead PDR, on the other hand, we find that the hy¬ 
drocarbon gas-phase chemistry is initiated by reactions 
of C + with CH leading to • Further reactions with 
H 2 lead to the formation of C 2 H + , C 2 H 2 " and 02 ^. 
The last two species recombine with electrons to form 
C 2 H and C 2 H 2 , respectively. C 2 H 2 can react with C + to 
form C 3 H+, but in fact the dominant formation route of 
CsH + at the edge of the cloud involves reactions between 
C 2 H and C + leading to C 3 ", followed by reactions with 
H 2 in the model. Once C 3 H + is produced, it reacts with 
H 2 to form the C 3 H 2 " and C 3 H 3 " ions which then recom¬ 
bine with electrons to form C 3 H and C 3 H 2 . CsH + can 
also recombine with electrons to form C 2 H, although the 
dominant formation route for C 2 H is the recombination 
of C 2 H 2 " with electrons: 


c 2 : 


CX 

A + 

-CH 



h 2 

c 2 h+ 

e — 

C2H2 

c 3 h+ 

Xc 3 H 2 






C+ H 2 .rf 


c 2 h+ 




C 3 H+ c 3 
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c 2 h- 

c + 3 


c 3 h+ 
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The dotted arrow in the scheme above marks a radia¬ 
tive association reaction, and the gray arrows indicate 
reactions with a temperature dependence of the rates. 
The CsH + hydrocarbon ion is thus a key chemical pre- 
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Spatially resolved /-CsH + emission in the Horsehead PDR 
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Fig. 3. — Photochemical model of the Horsehead PDR. (a) Den¬ 
sity profile nn = n(H) + 2n(H2). (b) Predicted abundance of small 
hydrocarbons, (c) Predicted abundance ratio relative to Z-CsH + . 
(d) Predicted intensity profile. Two lines are included for C 2 H: 
N =1-0 ( E u /k -4 K) and N = 3 - 2 (E u /k ~ 25 K); and for 
c-C 3 H 2 : J = 2-1 (E u /k ~ 6 K) and J = 3 — 2 (E u /k ~ 16 K). 
The observed intensity profile is overlaid with error bars. 


cursor of the small hydrocarbons in the gas phase. In 
particular, one would expect CsH + and C3H2 to have a 
similar spatial distribution if gas-phase chemistry alone 
is responsible for the observed C3H2 abundance. 

Fig. [ 3 ] displays the results of the photochemical model 
convolved to a resolution of 6" to facilitate the compar¬ 
ison with observations. The steep density profile of the 
Horsehead is shown in the upper panel (a). The abun¬ 
dances with respect to total hydrogen nuclei are shown 
in panel (b). The abundance ratios with respect to l- 
C 3 H + are shown in the panel (c). The observations are 
shown with error bars. In general, the abundances are 
well-reproduced by the pure gas-phase model at the PDR 
position (Ay ~ 1 . 5 ). At the PAH position (Ay ~ 0 . 05 ), 
on the other hand, the C2H and 0C3H2 abundances are 
underpredicted by an order of magnitude. The same con¬ 
clusion can be drawn from the abundance ratios, where 
the match between observations and model is better at 
the PDR than at the PAH position. We note that at the 
PAH position, the model better reproduces the observed 
/-C 3 H+ abundance than that of C2H and C3H2. This 
suggests the need of an additional formation mechanism 
to explain the observed hydrocarbon abundances at the 
edge of the cloud. 

The bottom panel (d) in Fig. [ 3 ] shows the modeled line 
intensity profile for each species and for different rota¬ 


tional transitions, normalized to their respective emission 
peaks. The observed profiles are overlaid with error bars. 
The modeled line intensity profiles of C2H and 0C3H2 
were computed by including the PDR model outputs (H2 
density, T^ n and column density) into the non-LTE ra¬ 
diative transfer code RADEX. In order to investigate ex¬ 
citation effects on the position of the emission peak, we 
included two transitions with different upper level ener¬ 
gies for each species: the 0C3H2 2i 2 — lo 1 (E u = 6 K) 
and C 2 H N = 1-0, J = 3/2 — 1 / 2 , F = 2 - 1 (E u = 4 K) 
(solid lines), and the 0C3H2 3i 2 — 3o 3 (E u = 16 K) and 
C 2 H N = 3 - 2 , J = 7/2 - 5/2, F = 3 - 2 (E u = 25 K) 
(dashed lines). These lines were chosen because their 
upper level energies are closer to that of the Z-C3H+ 
J = 5 — 4 line ( E u /k ~ 16 K). Because there are no 
available collisional coefficients for l-CsR + , we assumed a 
constant excitation temperature along the cloud to com¬ 
pute its line intensity profile. In the model all hydrocar¬ 
bons peak at the same position (5x ~ 15"). However, 
our observations show that the /-C3H + emission peak is 
shifted compared to that of the other hydrocarbons. We 
note that the emission peak of C2H and 0C3H2 does not 
change significantly between the lower and higher energy 
transitions. We have checked that including a constant 
excitation temperature for C2H and 0C3H gives simi¬ 
lar results. This suggests that the excitation tempera¬ 
ture gradient in the Horsehead takes place at very small 
spatial scales, and that the observed shift in the Z-C3H+ 
emission peak compared to the other hydrocarbons is not 
due to an excitation effect but reveals a real difference 
in the spatial distribution of their column densities that 
is not predicted by the current pure gas-phase chemical 
models. 

In general, the observed integrated intensity profile of 
all species is much broader than what the model predicts. 
In particular, the model underpredicts the Z-C 3 H+ emis¬ 
sion in the first PDR layers (5x < 10"). We run models 
with constant densities (10 4 — 10 5 cm -3 ) to explore the 
effect of the density profile on the /-C 3 H + intensity in 
the surface layers. We find that a higher density would 
slightly improve the agreement between model and obser¬ 
vations, but it would not reproduce the observations of 
many other tracers (e.g., HCO, H 2 , and dust continuum 
emission). Therefore, an additional formation mecha¬ 
nisms such as PAH photodestruction contributes to the 
formation of Z-C 3 H+ in these surface layers. Indeed, it 
has been suggested that C 2 H 2 , which would enhance the 
C 3 H+ abundance through reactions wi th C + , cou ld be 
a pro duct of photochemistry of PAHs (|Bierbaum et al.l 

Ilm!). 

4. CONCLUSIONS 

We have presented high-angular resolution ( 6 ") obser¬ 
vations of the hydrocarbon ion Z-C 3 H+ in the Horsehead 
photo-dissociation region. The Z-C 3 H+ emission is con¬ 
centrated toward the surface edge of the nebula, close 
to where other hydrocarbon chains, such as C 2 H and c- 
C3H2, and the 7.7 /i m PAH emission peak. However, in 
contrast to C 2 H and c- C3H2, Z-C 3 H+ is only abundant 
in the first PDR layers with almost no emission deeper 
inside the cloud. /-C 3 H + is thus a good tracer of hydro¬ 
carbon photo-chemistry. Moreover, /-C 3 H + peaks ^ 5" 
further outside the PDR than C 2 H and 0C3H2, which 
cannot be explained by an excitation effect. This is in 

































6 


V.V. Guzman et al. 


contrast to what current gas-phase chemical models pre¬ 
dict for the Horsehead physical conditions, that is, that 
all the hydrocarbons should peak approximately at the 
same position as their chemistry is closely linked. 

The inferred /-CsH + abundance at the PDR from the 
new high-r esolution o b serva tions is consistent with that 
derived bv iPetv et al.l (2012). In addition, we have com¬ 
puted the abundances of C2H, C3H2 and /-C3H + at three 
characteristic positions in the cloud. We find that a 
gas-phase chemical model can reproduce the hydrocar¬ 
bon abundances, including the /-CsH + abundance, at the 
PDR but it underpredicts the C2H, c- C3H2 and /-CsH + 
abundances at the PAH position, z.e., at the edge of the 
cloud or the first PDR layers. The fact that the disagree¬ 
ment between model and observations occurs at the edge 
of the cloud, where the PAHs emission peaks, is con¬ 
sistent with a top-down hydrocarbon chemistry, where 
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